Can the exciton— polariton regime be defined by its quantum properties? 
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Using a simple fully quantum model in an effective exciton scheme that takes into account the 
system-environment interaction, we study the different regimes arising in a microcavity-quantum 
dot system. Our numerical calculations of the emission linewidth, emission energy, integrated in- 
tensity and second- and third-order correlation functions are in good qualitative agreement with 
reported experimental results. We show that the transition from the polariton-laser to the photon- 
laser regime can be defined through the critical points of both the negativity and the linear entropy 
of the steady state. 
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The solid-state realization of Bose-Einstein Conden- 
sation (BEG) has been achieved in exciton-polariton 
systems 0, HJ. These quasi-bosons, arising from the 
strong coupling between photons and electron-hole pairs 
in semiconductor microcavities (/x-C), have a high criti- 
cal temperature due to their small effective mass (eight 
orders of magnitude smaller than hydrogen atom mass). 
After more than two decades of theoretical and exper- 
imental investigations, nowadays it is understood that 
the large occupation of the polariton ground state can- 
not be identified with usual thermodynamic equilibrium 
BEC states [H, 0, |] . Instead, the corresponding experi- 
mentally observed regime has been called polariton laser 
0, H, E| because of its dynamical nature and the gain 
in the light-emission intensity. The transition from this 
regime to a second one, identified with the well-known 
photon laser, has also been observed 

Cavity polariton systems have been studied from dif- 
ferent theoretical perspectives. Assuming thermal equi- 
librium, a trial wave function that takes into account the 
coherence properties of both light and matter, has been 
able to predict multiple phase transitions d, Q . On the 
other hand, when the matter-light state is obtained from 
some equation of motion (mean field dynamics lol- 12\, 
master equation in multiexcitonic scheme [l3[ , dissipative 
Jaynes-Cummings model (l3-21|). the dynamical char- 



acter of the polariton laser regime is conspicuous, and 
the non-Gibbsian character of the stationary state of the 
system is revealed. However, the elucidation of the mech- 
anisms behind the appearance of the different observed 
regimes is still an open problem. Our work is a first step 
in this direction. 

The aim of this Letter, the identification of the regimes 
observed in current experiments from the quantum prop- 
erties of the steady state, is possible due to the simplic- 
ity of our model. Indeed, we consider a single pumped 
radiator interacting with a leaky mode of the electro- 
magnetic field, ignoring collective effects. Our calcula- 
tions, however, are in qualitative good agreement with 



the experimental results. Moreover, we can correlate the 
entanglement, mixedness and the coherence functions of 
the steady state not only with the observed regimes but 
also with important physical parameters like the pump- 
ing rate and the detuning. In addition, we are able to 
provide a criterion to identify the "best" polariton that 
can be sustained by the system. 

We model a quantum dot (QD) embedded in a p-C, as 
a two- level system (ground \G) and excited \X) states). 
Its interaction with a single electromagnetic mode of fre- 
quency we, in the dipole and rotating wave approxima- 
tions, is described by the Hamiltonian (K = 1): 



H = idea' a + (ujc — A)<t^<t + g(aa^ 



a) a) 



(1) 



The detuning A is the difference between cavity mode 
and exciton energies, g is the matter-light coupling con- 
stant, a = \G) (X\ is the QD ladder operator, and a? 
(a) is the usual creation (annihilation) operator of the 
cavity mode. The Hamiltonian ([1]) commutes with the 
excitation number N = N p t + N ex — a^a + era] hence, 
it only causes transitions between matter-light states of 
the same excitation manifold. Polaritons are defined as 
the energy eigenstates H, and are explicitly given by 



\n,+) = sin$„ \G,n) 
In,-) = cos$„ \G,n) 



cos $„ | A, n— 1) 
-sin$ n |A,n-l), (2) 



where {\n)} denotes the Fock number states of the 
field and tan2$„ = 2<?^/n/A. We include two non- 
conservative processes, the loss of photons in the p-C 
(k) and the continuous pumping of excitons (P), in the 
master equation for the density operator p of the system 



^ = i[p, H] + \P(2^pa - aa^p - paa^) 
at z 

+ ^K{2apa) — a)ap — pa) a), 



(3) 



where we have made the Born-Markov approximation. 
We neglect other system-environment interaction mech- 
anisms (e.g., spontaneous emission, dephasing, photon 



2 



pumping, polariton pumping, etc.) because their effect 
is either small or is already effectively contained in the 
master equation. If a better adjustment with the exper- 
imental results is desired those mechanisms can be in- 
cluded and fitted [21| , but the qualitative physical image 
remains essentially unchanged. 

The basic assumption behind our approach, which fo- 
cuses on the steady state p ss of the equation of motion 
([3|), is that polariton lifetime is much longer than the 
time required to reach the asymptotic solution Q. The 
steady state p ss = p(n, P, g, A) of the system, is a func- 
tion of the dissipative rates k and P, the matter-light 
coupling constant g and the detuning A. In the weak- 
coupling regime g <C P, k, the pumping keeps the QD in 
its excited state while the dissipation steers the electro- 
magnetic field to its ground state. Other states are not 
significantly populated because matter excitation cannot 
be converted into photons. On the other extreme, the 
ultra-strong coupling regime g 3S> P, k, A, the long-time 
density matrix becomes (almost) diagonal in the basis 
of bare states \G/X,n). The larger the coupling g the 
smaller the difference of the populations of \G, n) and 
\X,n—l) (oc 1/ff 2 )- The coherences \ponXn-i\, which 
decay as 1/g, also vanish as the coupling g increases. 

In this work we focus on the strong-coupling regime 
g 3> P, k, in which the coherences PGnXn-i are small, 
but different from zero. In resonance they are purely 
imaginary. For small detunings they acquire a small real 
part. If the detuning increases, A> j, the matter-light 
interaction becomes dispersive, i.e., the energies of the 
matter states depend on the number of photons, and 
the mechanism which converts matter excitations into 
photons is suppressed. Thereby, large detunings corre- 
spond to a weak coupling regime. We conclude that in 
the regime \A\ ^ g ^> P, k, the steady state of the sys- 
tem is expected to exhibit a polaritonic behavior. Unless 
stated otherwise, the steady state solution p ss of |3]) is 
obtained for the initial condition p(0) = \G0) (G0\, and 
setting ujc — 1 eV, g = 1 meV and k = 5 x 10~ 2 meV, 
while A and P are varied in ranges similar to those of 
current experiments [!, 0, H3] • 

Evidence of the spontaneous coherence build-up as- 
sociated with polariton states are currently detected 
through the photoluminescence properties in quantum 
wells (QW) [Ei,@. We compare our theoretical pre- 
dictions in QDs with the experimental findings in QWs 
because: i) due to experimental difficulties no analogous 
results for QDs have been reported and ii) it is reasonable 
to assume that some of the physical mechanisms behind 
the exciton-polariton laser regime are the same in both 
cases. Additionally, the present approach may shed light 
on the separation of the collective effects in QWs from 
those of a QD single emitter. 

The calculated emission lincwidth, emission energy, in- 
tegrated intensity and number of photons are shown in 
fig. [T]as a function of the pumping rate. This numerical 
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FIG. 1. (color online), (a) Emission linewidth, (b) emission 
energy, (c) integrated intensity (continuous line) and aver- 
age number of photons (dashed line) and (d) second- and 
third-order correlation functions versus the incoherent exci- 
ton pumping P, for k — 5 x 10 -2 meV and A = 2.5 meV. 
The marked regions correspond to the polariton-laser and 
photon-laser regimes. 



calculation used the quantum regression theorem 
llq . [I8[ and the integrated intensity corresponds to the 
area under the curve of the peak associated to the tran- 
sition between the polaritons in the manifolds n — 1 and 
n = 2. The linewidth (fig HI a) exhibits the character- 
istic reduction in the polariton regime, and the subse- 
quent growth and decrease [H, 0, @ . The emission energy 
blueshifts (figUb) from the exciton transition frequency 
to the cavity mode frequency. In the polaritonic region 
the blueshift is smaller than observed in QWs 0, @] . In 
the intermediate region, f0 _1 meV < P < f meV, the 
calculated blueshift grows faster than the measured one. 
Despite the small slope changes in the integrated inten- 
sity as a function of P (figfflc) our results are consistent 
with the previous prediction of absence of threshold in 



a one-atom laser 23[. However, the nonlinearity of the 



model is evident in the curve of the average number of 
photons. Only for large values of the pumping (P > 1 
meV) , this curve is parallel to that of the integrated emis- 
sion intensity (whose slope is approximately one). 

Statistics of the emitted light can be characterized 
by the normalized second- g^(0) = (a^a^aa) / {a)a) 2 
and third-order g^(0) = (a^a^a^aaa) / (a^a) 3 coher- 
ence functions, plotted in fig.([TJd). For small pumping 
power (P < 2 x 1(T 2 meV) S (2) (0),5 (3) (0) < 1, a foot- 
print of quantum-like light, i.e. the partial state of the 
field is nearly a Fock state with small number of pho- 
tons, as expected. Indeed, our calculations show that in 
this regime p ss w \G) (G\ ® {sin<^ |0) (0| + cos^|l) (f |}, 
with ip sa 0. For intermediate values of the pumping 
(1CP 2 meV < P < 1CP 1 meV) both correlations func- 
tions monotonically grow beyond one (inset fig.l.d), as 
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FIG. 2. (color online). J\f(p as ) (continuous red line), Sl(Pss) 
(dashed black line) and v (dashed-dotted blue line) as a func- 
tion of the incoherent exciton pumping P, for k — 5 x 10~ 2 
meV and (a) A = 0, (b) A = 1 meV, (c) A = 2 meV, (d) 
A = 3 meV, (e) A = 7 meV and (f) A = 10 meV. Note that 
for A / the maximum of Af(p BS ), the minimum of Sl{Pss) 
and the inflection point of v coincide. 



has been experimentally observed [24(. We analyze this 
behavior below in the text. For larger pumping rates 
(P > 2 meV), the state of the field becomes coherent up 
to third order (when g^iO) — 9^(0) = 1), and thus the 
linewidth falls (fig.l.a). In this region, the state of the 
field obtained by the partial trace over the excitonic de- 
grees of freedom, has a fidelity of more than 0.99 with the 



random-phase coherent state |£ |ae l< ^) (ae^|, where 
\ae 1 ^) is an usual coherent state with \a\ 2 average num- 
ber of photons. This type of state has been proposed to 
describe the features of a (true) photon- laser [251 . 

We stress that the regimes identified through the char- 
acteristics of the emitted light, mirror quantum prop- 
erties of the system state. We describe those quantum 
properties by entanglement and mixedness of the steady 
state. Linear entropy and negativity are employed to 
quantify mixedness and matter-light entanglement, re- 
spectively. The former, defined as Sl(p) = 1 — Tr/> 2 , 
vanishes for pure states and is maximum for maximally 
mixed states. The latter is defined as Af(p) = 2 X)a<o \M> 
where A denotes the eigenvalues of the partial transpose 
of p [2^, 27 1 . Finally, searching for a relation between 
the energy of the system and its quantum characteris- 
tics (such as entanglement), we introduce the differential 
energy per excitation, v(k, P) — \d (H) jd (N), as a con- 
venient measure of energy per particle. The factor of 2 
in the definition of v has been chosen to satisfy the con- 
dition v(P <C g, K <C A) w bJc — A. The negativity, 
linear entropy and the differential energy per excitation 
are depicted in fig. [2] 

Assuming the strong-coupling regime, the formation 
of the polariton is hindered by two different mechanisms, 
which depend on the detuning. While the mixedness of 



the state is large for small detuning (fig. [2]a and[2jb), 
matter and light decouple for large detuning (fig. [2jf). 
In the intermediate region g < |A| < lOg, where the po- 
lariton is well defined, both mechanisms compete. The 
matter-light entanglement is enhanced, and the entropy 
of the asymptotic state of the system decreases with in- 
creasing detuning. The negativity of the steady state 
Af(f>ss) vanishes for small and large values of the pump- 
ing power, and attains a maximum at the point P = n 
-the mid-point of the polariton-laser region. 

Now, we are able to give a possible explanation of 
the behavior of the correlation functions in the po- 
lariton regime. The asymptotic state in the polari- 
tonic region is a mixed entangled state which satis- 
fies 3^ 3) (0) > 3^ 2) (0) > 1. This behavior of the cor- 
relation functions is a rather generic feature. As an 
example, we consider another mixed entangled state 
P P oi(n) = J2n Pn{n) \n, +) (n, +|. Since the probabili- 
ties P n (n) = e~ n n n / (n\) are Poisson weights, this is a 
polariton coherent state. However, the reduced photon 
state has super-poissonian statistics, i.e., the second- and 
third-order coherence functions can not be expected to 
be unity. Hence, matter-light entanglement is a viable al- 
ternative to the standard explanation of the unexpected 
behavior of these functions, based on polariton-polariton 
and polariton-phonon interactions. 

The differential energy per excitation v is plotted in 
fig. [2] as a function of P for several values of the de- 
tuning. For A = 0, v is independent of P and equals 
to the cavity mode energy, since the number of photons 
becomes much larger than the number of matter excita- 
tions. Our numerical results show that for |A| > g, and 
small (P < 1(T 2 meV) or large (10 meV > P > 10" 1 
meV) pumping rates, v is almost a constant, equal to the 
exciton (photon) energy in the former (latter) case. The 
same constant values would had been obtained with the 
definition v = H/N. For intermediate values of P, v dis- 
plays an inflexion point at P — k, where v = loc — (A/2) 
is halfway between the exciton and photon energies. 
Moreover, since the mean number of excitations is one, 
it is tempting to define P — n as the condition for the 
"optimum" polariton. In order to quantify this idea we 
compare the steady state with the polariton states |n, ±) 
defined by ([2|), using the sequence of non-zero fidelities 
F n ± = (n, ±\p ss \n, ±). For small values of P (figOa), 
when pgogo is much larger than the other populations, 
only F\- does not vanish -however it is relatively small-. 
For P = k (figElb) F]_ increases up to more than 0.95, 
while the remaining fidelities are still small. Hence, the 
steady state p ss is quite similar to the Ai-polariton |1— ). 
As P increases, F n ± is non-zero for larger excitation- 
numbers, but their values are very small. 

The excitation number (N) symmetry associated with 
the Hamiltonian (JlJ is broken in the time evolution pro- 
vided by the master equation ([3]), in the sense that the 
asymptotic state of the system cannot be labeled with a 
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FIG. 3. (color online). Sequence of non-zero fideli- 
ties F n ± between the steady state p S3 and the A n - 
lower(black)/upper(green) polaritons \n, ±) for k — 5 x 10~ 2 
meV, A = 3 meV and (a) P = 2 x f 0" 3 meV, (b) P = 5 x I(T 2 
meV, (c) P = 1 meV and (d) P = 5 meV. |A n j denotes the 
excitation number of the polariton manifold A n . 

single eigenvalue of N. If polariton-like behavior is actu- 
ally present, restoration of the symmetry is expected. In 
order to account for this effect, we introduce the partici- 
pation ratio PR = J2^=o ^n> where- P n is the probability 
to have n excitations in the asymptotic state. This quan- 
tity, which varies from zero -all excitation numbers are 
cquiprobable- to one -only one occupied manifold-, dis- 
plays a global maximum at zero pumping rate and a local 
maximum at P ~ k, in the strong coupling regime, sig- 
naling a partial restoration of symmetry. This can be 
understood as a combined effect of the decrease of the 
mixedness of the state and the increase of its entangle- 
ment, occurring at P w discussed above (again, 
A~ 5 ). 

As we have seen, the optimum polariton exhibits max- 
imum negativity, minimum linear entropy, a local maxi- 
mum of the participation ratio and an inflection point of 
the differential energy per excitation, provided that the 
system is in strong coupling. To quantify our previous 
qualitative arguments -which show that polaritons can- 
not be sustained neither for small nor for large detunings- 
, it is worth examining the behavior of the linear entropy 
and the negativity at P = function of the detun- 

ing (figH]). Three regions can be identified. In the first 
region, | A| < 0.68 meV, the negativity of the steady state 
is exactly zero and the linear entropy is larger than 0.72. 
In the second region, 0.68 meV < |A| < 10 meV, the lin- 
ear entropy and the negativity are still significative. In 
the third region the dissipative polariton is very close to 
the Hamiltonian polariton |1,±), for the corresponding 
A. Nevertheless, the dressed states ^ are virtually the 
bare states. 

Our results might provide a guide to experiments, in 
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FIG. 4. (color online). Linear entropy (dashed black line) 
and negativity (continuous red line) of the steady state of 
the system p as , and negativity of polaritons of the excitation 
manifold Ai (dashed-dotted blue line) as a function of A for 
P = /t = 5x 10" 2 meV. 

the sense that, in the strong-coupling regime, both the 
pumping rate and the detuning have to be carefully ad- 
justed. In our model, the best matter-light correlation 
properties occur at P = n and A 3g, where N ~ 0.32 
and Sl ~ 0.29. From the theoretical point of view we 
propose the following criterion: if the negativity rises 
above 0.25, close to its maximum possible value, its in- 
flection points, as a function of the pumping power, can 
be used to define the polaritonic regime. When this con- 
dition is fulfilled all the quantifiers that we have exam- 
ined exhibit a characteristic change. The emission energy 
presents a blueshift, the differential energy per excitation 
has an inflection point, the emission line decreases, the 
second- and third-order correlation functions increase be- 
yond their value for photon coherent states, the entropy 
decreases and the negativity increases. With the excep- 
tion of the first two, these changes can be understood as 
a coherence gain of the asymptotic state of the system. 
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